Thermal Casimir vs Casimir— Polder forces: Equilibrium and non-equilibrium forces 
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We critically discuss whether and under what conditions Lifshitz theory may be used to de- 
scribe thermal Casimir-Polder forces on atoms or molecules. An exact treatment of the atom-field 
coupling reveals that for a ground-state atom (molecule), terms associated with virtual-photon ab- 
sorption lead to a deviation from the traditional Lifshitz result; they are identified as a signature 
of non-equilibrium dynamics. Even the equilibrium force on a thermalized atom (molecule) may be 
overestimated when using the ground-state polarizability instead of its thermal counterpart. 

PACS numbers: 34.35.+a, 12.20.-m, 42.50.Ct, 42.50.Nn 



Dispersion forces such as Casimir and Casimir-Polder 
(CP) forces are of increasing relevance in nanophysics 
[l| ; recent successes include the chemical identification of 
surface atoms via atomic force microscopy Q as well as 
the construction of novel biomimetic dry adhesives [3|. 
They become important in efforts to miniaturize atom 
chips [i| and play a key role in experiments placing upper 
bounds on non-standard gravitational forces [5[ . 

In all of these thorough understanding of dis- 

persion forces under realistic conditions must account for 
their temperature dependence induced by thermal pho- 
tons p|. A series of high-precision experiments [3] has 
triggered a renewed interest in this thermal Casimir force 
[H by opening the perspective of its experimental inves- 
tigation [9]. It was noticed that, depending on the model 
chosen to describe the metal response, Lifshitz theory can 
yield different answers for the temperature dependence of 
the Casimir energy between two metal plates |lo| . The 
resulting debate concerning the correct description of the 
thermal Casimir force [ll| will ultimately have to be set- 
tled by experiments. Non-equilibrium situations of two 
plates of different temperatures have recently been sug- 
gested as sensitive probes to the quantum electrodynam- 
ics (QED) of the Casimir effect 

The CP force on single atoms can be measured indi- 
rectly via spectroscopic means; clear evidence for ther- 
mal frequency shifts has been found for atoms inside a 
cavity [l3| and their signature has been detected in the 
interaction of atoms with a sapphire surface [13 ]. The 
novel non-equilibrium CP forces predicted for the case 
of different surface and environment temperatures [IH 
have recently been observed via their effect on the center- 
of-mass oscillations of a trapped Bose-Einstein conden- 
sate [ifj. While some theoretical studies of the thermal 
CP force are based on a QED treatment of the atom- 
field coupling 17, 3], the vast majority of investigations 
invokes a macroscopic calculation using Lifshitz theory 

(j [Til |20| or a linear-response description of the atom 

15l 2lT | (leading to equivalent results) . 

The macroscopic approach to the CP force is based on 
a very close relation between Casimir and CP forces. It 



is the validity of this one-to-one correspondence and its 
results for the CP force that we intend to investigate in 
this Letter, so let us briefly recount the argument in its 
traditional form: Generalizing the famous Lifshitz result 
for two plates [f|, the thermal Casimir force on a homo- 
geneous body of arbitrary shape occupying a volume V 
in free space due to the presence of another body can be 
given as a Matsubara sum [22J, l23j 

r. OO 

F = 2k B T / d 3 r V (1 - ±6 N0 ) 

Trjl x [V x V x +(Wc) 2 ] G« (r, r', i&v) x V') 



(Wc) 2 V-G (1 >(r,r,i6v) 



(1) 



with Matsubara frequencies £jy = 2irk B TN /h. Here, G*- 1 -* 
is the scattering part of the classical Green tensor 

[V x V x - (w/c) 2 e(r, w)] G(r, r', w) = 5{r - r') (2) 

associated with the bodies that are characterized by their 
dielectric permittivity e(r,oj). The result ([1]) can be ob- 
tained from the zero-temperature force by applying the 
simple replacement rule 



oo 

2k B T (! 

N=0 



IM/^) (3) 



which is equivalent to replacing the zero-point energy 
by a thermal spectrum [n(cj) + i]fiw with pho- 
ton number n(uj) = I / [e Tmj I - l] . To derive the 
CP force, one assumes that one body consists of a di- 
lute gas (number density rf) of atoms with polarizability 
a{oj), so that its permittivity may be approximated by 
e{lo) = 1 + r/a(w)/eo. After a linear expansion in a, it 
follows that the thermal CP force is given by 22, 23| 



F(r A ) = -fi k B T IM&^&v) 



N=0 



Y4TrG (1) (r4,r-A,i6v). (4) 
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As an immediate consequence of its macroscopic deriva- 
tion, this force expression is very similar in structure to 
the respective Casimir force. In particular, the replace- 
ment rule © is again valid which often serves as a start- 
ing point for studies of the thermal CP force 24( . 

In view of the high precision of modern experiments, 
the current widespread use of the macroscopic result (|4|) 
and the replacement rule ((3]) for calculating thermal CP 
forces calls for a critical discussion of their validity. In 
this Letter we clarify the extent to which Eq. ((4| can 
be used on the basis of a direct calculation of the time- 
dependent thermal CP force from the atom-field cou- 
pling Hamiltonian. We will demonstrate that force com- 
ponents associated with thermal-photon absorption will 
generically lead to deviations from Eq. (j4|). As we will 
show, the macroscopic result provides a reasonable ap- 
proximation to the force on a fully thermalized atom, 
provided that the atomic polarizability is correctly inter- 
preted. 

In electric dipole approximation, the force on an atom 
or a molecule (called atomic system in the following) pre- 
pared in an incoherent superposition of internal eigen- 
states |n) is given by [23I |2{| 



F(r A ,t) = ([Vd-E(r)] >, 



(5) 



where it is assumed that the internal and translational 
motion of the atomic system decouple in the spirit of 
Born-Oppcnhcimcr approximation. In order to evaluate 
this expression, one needs to solve the coupled dynamics 
of the electromagnetic field and the atomic system, which 
is governed by the Hamiltonian [26| 

//•oo 
d 3 r J dujhujf\r,uj)-f(r,uj) 

+ ^E n A nn -^A mn d mn -E(r A ) (6) 



(d mn — (m\d\n), A mn = \m){n\). The bosonic dynami- 
cal variables and /, which describe the elementary ex- 
citations of the body-assisted electromagnetic field, can 
be used to construct an expansion of the electric-field 
operator E(r) = f^°dujll(r, u) + H. c. according to 



E(r, w) = i y/h/^eo) (u/c) 2 J dV ^Jme(r',oj) 

xG(r,rV)-/»). (7) 

The expansion coefficients are given in terms of the clas- 
sical Green tensor ([2]) which obeys the integral relation 

Ll> 2 f 

— / d 3 s Ime(s,(j)G(r,s,w)-G*(s,r',w) 

= ImG(r,r». (8) 



The Hamiltonian ([6]) implies that the atom-field dy- 
namics is governed by the system of equations 



7- ^2(d nk A mk - d km A kn )-E(r A ), (9) 



f(r,u) = -iwf(r,u) + y / Ime(r,uj)/(hire ) (w/c) 2 

x A mn d mn -G*(r A ,r,u). (10) 



We eliminate the field by formally integrating Eq. (flT 
and substituting the result back into Eq. ([9]), which we 
arrange in normal ordering. After using the integral re- 
lation ([5]), we obtain 

^mn(^) i^mra^-mn (^) ~t~ ^mn(^) 
i f°° 

+ &u{e lwt E\r A ,u)){d nk A mk {t) - d km A kn (t)] 

k ^° 

+ eT m [A mk (t)d nk - A kn (t)d km ] -E(r A , w)} (11) 
where 

Z mn {t) — ^ / dujuj 2 d k rlmG(r A ,r A ,uj) / dr 
j,k,l Jo Jt) 
.{e-^-^ [d jm A jn (t) - d n] A mj (t)]A kl (T) 

+ e^-^Auir) [d nj A m3 {t) - d jm A jn (t)] } (12) 

denotes the zero-point contribution to the internal atomic 
dynamics. The thermal contribution can be determined 
iteratively by substituting the self-consistent solution 



x [A mk (T)d nk - A kn (r)d km ]-[E(r A ,co)e iajr + H.c] 

(13) 

into the truncated Eq. (jlip [without the zero-point con- 
tribution Z mn (t)] back into itself and taking expectation 
values. With the field initially being prepared in a ther- 
mal state p T = exp[- J ff F /(fc B T)]/Tr{exp[- J ff F /(fc B T)]}, 
Hf = J d 3 r duj hjjy (r,w)-/ (r,w), the non-vanishing 
averages of the field operators are 

= {hpL Q /iT)n{uj)uj 2 Im G(r, r', u)6(u - u') (14) 

[recall Eqs. ([7]) and ©]. One thus obtains a closed system 
of equations for the internal atomic dynamics 

{A mn j= iUmn{A mn ^ + (Z mn ^+(T mn ) , (15) 
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with the thermal contribution being given by 

(■^nn(^)) 



k,l 



x {e ua »»( t - r 5d nfc -ImG(rA,rA,w) 
• [d irn (A ;fe (T))-d fc; (A m/ (T))] 
+ e i ^" (t - T) d fem .ImG(rA,rA,w) 

•[^(ifeKr))-^^^^))]}. (16) 

For weak atom-field coupling, this system can be solved 
in Markov approximation using the relations \A mn {ry) ~ 

e-' lC>mn{t - T) (A mn (t)) and / f dre ia: ( f - T ) ~ 7n5(x) + iP/x. 

For a non-degenerate system, the off-diagonal elements 
of the atomic density matrix a decouple from one another 
as well as from the diagonal elements, and one finds that 
the internal atomic dynamics follows the rate equations 



(17) 



ow(t) = [— iti-mn - |(r m + T n )]a mn (t), m^n (18) 
(<7 m „ = (m\<j\n) = (A nm )). Here, the total loss rates read 

r„ = ^r„ fc = ^^2tf lk {e(Cj nk )[n(Cj nk ) + 1] 

k k 

+ Q(uj k n)n(oJkn)}d n k-lmG(r A ,r A , \Co nk \)-d kn (19) 



and the shifts of the atomic transition frequencies u> v 
Lo mn + 5uo m - 5uj n are given by 



Sw n = ^ Suj nk = — \^ V 

k k 

ti(lu) 



duj lj' 



n{u) + 1 
u) n h — ui 



Unk. + OJ 



dnk-IaxG(rA,rA,w)-dkn', (20) 



they reduce to the perturbative result derived in Ref. [16 
upon setting uj nk — w nk in the denominators. 



With these preparations, we evaluate the CP force ([5|) 
with the help of Eqs. 0, (13]) and (O and the solution 
to Eq. (JTOl) as 



F(r A ,t) 

= —^2 / dajLo 2 Vd nk -lmG {1 \r,r Al uj)-d kn \ r=j 



n . k 



x / dr(i„. n (r)){nHet i ^"''-( r " +r ^ 2 l( t ~ T ) 
Jo 

+ [n(w) + l]el- f (— « »*)-(r„+r fc )/a](t-r) j + c c (21) 



Here we have used the correlation function 



{A mn {t)Au{r)) = 5 n k(A m i(r))( 



„-(r m +r n )/2](t-r) 
(22) 

which follows from Eq. (| 18() by means of the quantum 
regression theorem [271 ]. Evaluating the r integral in 
Markov approximation and the lo integral by means of 
contour-integral techniques (cf. Ref. [25[), one finds that 
the thermal CP force on an atomic system prepared in 
an incoherent superposition of energy eigenstates is given 
by F(r A ,t) = Y, n v n n(t)F n (r A ) with force components 



F n {r A ) = -fi k B T C 1 - §M& VTr {K(i6v) + a n (-i^ N )]-G^(r A ,r,^ N )}\ r=rA + ^^{O^)^ 

JV=0 k 

x [n(n nk ) + l]Vd nk -G^(r,rA,n nk )-d kn \ r=rA -e(Q kn Mln(QUVd nk -G^\r,r Ai ntJ-d kn \ r=rA +c.c^ (23) 

I 

[f2 n fe = u)„fe + i(r„ + Ffc)/2] and atomic/molecular polar- lated in Ref. [25| . 
izability 

In order to compare this force with the macroscopic 
result ([31) , we consider an isotropic atomic system in the 
This result generalizes the zero-temperature force calcu- perturbative limit £l nk ~ ui nk , whereby Eq. (|23|) simpli- 



a n (uj) 



-y 



dnkdkn 



dkndnk 



-Kk + u 
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fi.es to 



_'JN0) 

x ^ N a n (^ N )V A TTG (1) (r A ,r A ,^ N ) 

+ y^nfc{6(^nfc)N^nfc) + 1] - 9(w fc „)n(u; fe „)} 

fc 

x |d„ fc | 2 V A TrReG (1) (r A ,r A , W „ fe ) (25) 

and coincides with the negative gradient of the frequency 
shift (gnj) given in Ref. [18]. In addition to the Matsub- 
ara sum, the thermal CP force has resonant contributions 
proportional to n(uj n k) and respectively, which 

are due to the absorption and emission of photons by 
the atomic system. Even the force Fq(t a ) on a ground- 
state atom or molecule exhibits resonant force compo- 
nents -^fi ujl n(uj k o)\do k \ 2 V A TiReG {1) {r A ,r A ,uj k o) as- 
sociated with thermal-photon absorption which lead to a 
discrepancy with the macroscopic result {4]). These res- 
onant forces can be observed on time scales which are 
short with respect to the inverse ground-state heating 
rates T^" 1 , and their magnitude scales with the number 
of thermal photons at the respective transition frequency 
n(wfeo)- Polar molecules thus present an ideal candidate 
for studying them, since their heating time can be of the 
order of several seconds, and the thermal photon number 
at the relevant vibrational and rotational transition fre- 
quencies in the microwave regime can reach values of up 
to a few hundred at room temperature 28[ . As discussed 
in Ref. the enhanced frequency shifts observed very 
recently via selective reflection spectroscopy of cesium in 
the far infrared might be due to a resonant thermal effect 
of this kind. 

Our approach allows us to discuss the full dynamics 
of the CP force. It is thus able to reveal that resonant 
force components associated with absorption and stim- 
ulated / spontaneous emission of photons are a genuine 
non-equilibrium effect: According to Eq. (fTTj) , the atomic 
system reaches thermal equilibrium with its environment 
in the long-time limit and is described by the thermal 
state a T = e-^^^/^^/Trfe-^n^^/CteT)] 

(E n = E n + h5uj n ). In thermal equilibrium, all reso- 
nant force components cancel and the force can be writ- 
ten in the form of Eq. (0| only if the atomic polariz- 
ability has been identified as the thermal polarizability, 
ar(w) = J2 n a T,nnC( n (uj). In the limit of a single dom- 
inant transition, this equilibrium force is smaller than 
that resulting from Eq. (j4j with ao(uj) by a factor tt 
= tanh[toi /(2fc B T)] = \2n{u) W ) + The force on 

ground-state atoms, which is dominated by electronic 
transitions, is insensitive to this effect at room tem- 
perature due to the very small thermal photon num- 
ber: For Rb (wio = 2.37xl0 15 Hz rad one obtains 
r T = 1 - 1.3 x 1(T 26 at T = 300K. As in the case of the 
non-equilibrium force discussed above, a noticeable devi- 



ation from the Lifshitz result is to be expected for polar 
molecules: For CaF, force components associated with 
vibrational transitions (luiq — 1.15 x 10 14 Hz rad [I(|) are 
reduced by a factor tt — 0.90, those associated with rota- 
tional transitions (u)\q = 1.32 x 10 11 Hz rad |3lJ) even by a 
factor tt = 0.0017 — they are thus strongly reduced with 
respect to the prediction from Lifshitz theory. 

In conclusion, a full quantum-mechanical treatment of 
the atom-field interaction has revealed that the thermal 
CP force on a ground-state atom or molecule cannot 
be obtained from a macroscopic calculation while the 
force on a fully thermalized atom can — provided that 
one uses its correct finite-temperature polarizability. As 
shown, the force on thermalized molecules is consider- 
ably smaller than suggested by Lifshitz theory. The res- 
onant contributions to the force that occur for a ground- 
state atom are due to the absorption of thermal photons 
and manifest a non-equilibrium effect. In particular, they 
may dominate the thermal CP force on polar molecules, 
where they present both a considerable limiting factor 
for the miniaturisation of molecular surface traps and a 
novel probe to the surface-assisted thermal dynamics of 
these molecules. The discussion presented here imme- 
diately generalizes to scenarios such as those suggested 
in Ref. [l5j where different parts of the environment are 
held at different temperatures. 

This work was supported by the Alexander von Hum- 
boldt Foundation and the UK Engineering and Physical 
Sciences Research Council. 



[1] 
[2] 
[3] 

[4] 

[5] 
[6] 



[9] 
[10] 

[11] 



H. G. Craighead, Science 290, 1532 (2000). 

Y. Sugimoto et. al, Nature 446, 64 (2007). 

M. T. Northen and K. L. Turner, Nanotechnology 16, 

1159 (2005). 

J. Fortagh and C. Zimmermann, Rev. Mod. Phys. 79, 
235 (2007). 

R. S. Decca et. al, Ann. Phys. 318, 37 (2005). 
E. M. Lifshitz, Sov. Phys. JETP 2, 73 (1956); I. E. 
Dzyaloshinskii, E. M. Lifshitz, and L. P. Pitaevskii, Avd. 
Phys. 10, 165 (1961). 

S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997); U. Mo- 
hideen and A. Roy, Phys. Rev. Lett. 81, 4549 (1998); 
R. S. Decca et. al, Phys. Rev. Lett. 91, 050402 (2003); 
for reviews, see M. Bordag, U. Mohideen, and V. M. 
Mostepanenko, Phys. Rep. 353, 1 (2001); S. K. Lam- 
oreaux, Rep. Prog. Phys. 68, 201 (2005). 
C. Genet, A. Lambrecht, and S. Reynaud, Phys. Rev. 
A 62, 012110 (2000); G. L. Klimchitskaya and V. M. 
Mostepanenko, Phys. Rev. A 63, 062108 (2001). 
M. Bostrom and B. E. Sernelius, Phys. Rev. Lett. 84, 
4757 (2000). 

M. Bordag et. al, Phys. Rev. Lett. 85, 503 (2000); J. S. 
H0ye et. al., Phys. Rev. E 67, 056116 (2003). 
V. M. Mostepanenko et. al, J. Phys. A: Math. Gen. 39, 
6589 (2006); I. Brevik, S. A. Ellingsen, and K. A. Milton, 
New J. Phys. 8, 236 (2006). 



■5 



[12] M. Antezza et. ai, Phys. Rev. Lett. 97, 223203 (2006); 

V. B. Bezerra et. ai, Eur. Phys. J. C. 52, 701 (2007). 
[13] M. Marrocco et. al., Phys. Rev. Lett. 81, 5784 (1998). 
[14] P. C. d. S. Segundo et. ai, Laser Phys. 17, 983 (2007). 
[15] M. Antezza, L. P. Pitaevskii, and S. Stringari, Phys. Rev. 

Lett. 95, 113202 (2005). 
[16] J. M. Obrecht et. al, Phys. Rev. Lett. 98, 063201 (2007). 
[17] T. Nakajima, P. Lambropoulos, and H. Walther, Phys. 

Rev. A 56, 5100 (1997); S.-T. Wu and C. Eberlein, Proc. 

R. Soc. Lond. Ser. A 456, 1931 (2000). 
[18] M.-P. Gorza and M. Ducloy, Eur. Phys. J. D 40, 343 

(2006). 

[19] M. Antezza, L. P. Pitaevskii, and S. Stringari, Phys. Rev. 
A 70, 053619 (2004). 

[20] J. F. Babb, G. L. Klimchitskaya, and V. M. Mostepa- 
nenko, Phys. Rev. A 70, 042901 (2004); E. V. Blagov, 
G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. 
Rev. B 71, 235401 (2005); G. L. Klimchitskaya, E. V. 
Blagov, and V. M. Mostepanenko, J. Phys. A: Math. 
Gen. 39, 6481 (2006); M. Bordag et. ai, Phys. Rev. B 74, 
205431 (2006); E. V. Blagov, G. L. Klimchitskaya, and 
V. M. Mostepanenko, Phys. Rev. B 75, 235413 (2007). 

[21] A. D. McLachlan, Proc. R. Soc. Lond. Ser. A 274, 80 
(1963); C. Henkel et. al, J. Opt. A: Pure Appl. Opt. 4, 



S109 (2002). 

[22] C. Raabe and D.-G. Welsch, Phys. Rev. A 73, 063822 

(2006) . 

[23] S. Y. Buhmann and D.-G. Welsch, Prog. Quantum Elec- 
tron. 31, 51 (2006). 

[24] J. N. Israelachvili, Q. Rev. Biophys. 6, 341 (1974); V. M. 
Nabutovskii, V. R. Belosludov, and A. M. Korotkikh, 
Sov. Phys. JETP 50, 352 (1979); M. Bostrom and B. E. 
Sernelius, Phys. Rev. A 61, 052703 (2000); T. H. Boyer, 
Phys. Rev. A 11, 1650 (1975). 

[25] S. Y. Buhmann et. al, Phys. Rev. A 70, 052117 (2004). 

[26] L. Knoll, S. Scheel, and D.-G. Welsch, in Coherence and 
Statistics of Photons and Atoms, edited by J. Pefina (Wi- 
ley, New York, 2001), p. 1. 

[27] W. Vogel and D.-G. Welsch, Quantum Optics (Wiley- 
VCH, Berlin, 2006), 3rd ed. 

[28] N. Vanhaecke and O. Dulieu, Mol. Phys. 105, 1723 

(2007) . 

[29] J. E. Sansonetti, J. Phys. Chem. Ref. Data 35, 301 
(2006). 

[30] M. Pelegrini et. al, Braz. J. Phys. 35, 950 (2005). 

[31] L. A. Kaledin et. al, J. Mol. Spectrosc. 197, 289 (1999). 



